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Abstract 

The attractors of Iterated Function Systems in Euclidean space - IFS fractals - have 
been the subject of great interest for their ability to visually model a wide range of 
natural phenomena. Indeed computer-generated plants are often modeled using 3D IFS 
fractals, and thus their extent in virtual space is a fundamental question, whether for 
collision detection or ray tracing. A great variety of algorithms exist in the literature 
for finding bounding circles, polygons, or rectangles for these sets, usually tackling 
the easier question in 2D first, as a basis for the 3D bounding problem. The exist- 
ing algorithms for finding bounding circles are mostly approximative, with significant 
computational and methodological complexity. We intend to hereby introduce explicit 
formulas for bounding circles in the plane, and some generalizations to space, thereby 
providing readily applicable bounding sets for IFS fractals.* 

MSC class: 28A80 (primary); 68U05, 52A27 (secondary) 
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Earlier efforts to provide bounding sets to IFS fractals, have been mostly algorithmic with 
significant computational complexity, while explicit bounding circles are yet to be intro- 
duced. Our circumcircles generalize the concept of the circumscribed circle of triangles in 
classical geometry to IFS fractals, hinting at a novel polygonal viewpoint of IFS fractals. 
Our general bounding circle resolves the bounding problem for any number of maps and 
dimensions. Our paper is intended for mathematicians and computer scientists alike. The 
latter may appreciate our discussion of the visualization of fractals and our ready bounding 
circle formulas, while the former, our attempt to take the bounding problem to a geometrical 
level. Indeed our primary goal is to reason the usefulness of bounding circles as a numerical 
panacea for any geometrical problem related to IFS fractals. 



*A preprint of the article has been submitted for consideration in the Journal of Computational Geometry. 
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The relevance of bounding sets lies in their property of "iterative containment" by which 
they provide improving approximations to the IFS fractal. The idea of "bounding" is im- 
portant so that no points of the fractal are excluded, and thus we may devise numerical 
theorems and algorithms for virtually any geometrical question about an IFS fractal, such as 
its intersection with Euclidean sets, or even other fractals. The initial bounding set can be a 
circle, the convex hull, or any other, which may then be iterated for better approximations 
to the fractal. 




Figure 1: A bounding circle iteratively approximating the fractal. 



The challenge is therefore to find bounding circles which are as tight as possible, and prefer- 
ably have a ready formula for their center and radius. In the current literature, such explicit 
bounding circles are scarce, while there are several defined via algorithms. Virtual plants 
as well as other objects in 3D, are often generated in an iterative way, essentially producing 
"embellished" IFS fractals. Their bounding can be relevant for multiple applied reasons. For 
morphing fractals, real-time bounding via explicit formulas is most ideal. Though research 
often begins by treating the 2D bounding problem first, the eventual goal is to generalize 
such results to 3D space. Our "general bounding circle" indeed also extends to 3D. 

For the visualization of IFS attractors, or really for any object in virtual space, there are two 
primary questions arising from the problem of lighting the object or "ray tracing" it: the ob- 
ject's spatial extent, and its intersection with various rays of light. With classical Euclidean 
objects these questions are much easier. As illustrated by Figure 2, ray tracing is composed 
of a series of steps for shading an object and its environment. Testing for whether a camera 
ray intersects the object involves: (1) sending out a test ray, (2) testing for intersection with 
the convex hull of the object, (3) if the ray intersects the hull, testing for intersection with 
the object itself. Determining the effect that a light source has on the object, involves similar 
steps. Paper [7] provides a full description of the rendering / visualization of 3D IFS frac- 
tals. It emphasizes the importance and difficulties of the ray / line intersection problem for a 
fractal. The paper also describes a bounding sphere, and utilizes it for the rendering process. 
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Other than ray tracing, the fundamental question of handling an IFS fractal in virtual 
space may arise. We must be able to calculate its interaction with other virtual objects. We 
must know its extent to be able to bump into it, or hold it in an avatar's hands even. Thus 
we may observe that the bounding problem leads to a resolution of fundamental geometrical 
questions in computer graphics applications, related to IFS fractals. 




Figure 2: Ray tracing. Figure 3: A ray-traced IFS tree. 



We will hereby review some of the main cornerstones in bounding circle research. Essen- 
tially, the main intentions here are to give some sort of approximation to a fractal, in order 
to assist its display. Article [16] ends with a bounding circle for two-map IFS, similar to 
our general bounding circle for multi-map IFS. That method however does not optimize for 
a center, but always uses the origin. Paper [4] comes even closer to our definition of the 
general bounding circle, having a variable center. However in this form, it is unclear how the 
optimal center may be found in general, possibly using an optimization algorithm. By our 
alternative formulation however, it becomes clear that the minimization involved is equiv- 
alent to the Smallest Circle / Ball Problem of finitely many points, for which specialized 
efficient algorithms exist. Paper [2] discusses the connectivity of a special class of two-map 
IFS fractals. Though not its primary focus, the paper introduces a bounding circle for a 
certain class of two-map IFS, which is in fact a special case of our circumcircle. In [13] we 
are introduced to a method for multi-map IFS, similar to our circumcircle method. In fact, 
the paper even conjectures the existence of the circumcircle, as an ideal tangential solution to 
the thereby presented numerical method, in the two-map case. Again this method gives no 
explicit formula, but assumes the use of an efficient optimization algorithm at each iterative 
step, leading to a bounding circle. The method presented in paper [3] gives an algorithm 
for finding not a single bounding circle, but rather a finite set of circles whose union covers 
the fractal. An outer bounding circle to this set of circles can thus be calculated and will 
be a bounding circle of the fractal itself. Considering the figures, the method seems rather 
inefficient, and the size of the resulting bounding circle is generally quite large relative to 
the fractal. Both [5] and [14] assume an initial good center estimate which remains fixed 
throughout the algorithm, and makes no improvement upon it. In [11] a novel algorithm is 
presented which is supposedly able to arrive at the absolute tightest bounding sphere of the 
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fractal, via a series of subroutines. The sphere however is restricted to be centered at the 
centroid. The method is founded upon the theory of IFS moments and the simplex method 
from linear optimization, that has exponential complexity in the worst case. Thus the imple- 
mentation of the method requires significant computational and methodological complexity. 
In our Remarks, we reason a much simpler method for finding the tightest bounding circle, 
by improving an existing one. 

We shall now turn to introducing the basic mathematical concepts related to IFS fractals 
and bounding circles. For an introduction to fractals, see [10, 1, 6]. 

1 IFS Fractals and Bounding Circles 

The attractors of Iterated Function Systems - IFS Fractals - are likely the most basic types 
of fractals possible, first defined in [8]. They are the attractors of a finite set of affine linear 
contractive maps - the "function system" - which when combined and iterated to infinity, 
converges to an attracting limit set, the IFS fractal itself. There is a fundamental motiva- 
tion in studying linear contractive maps in the plane, as they are the simplest contractions 
possible, and give rise to many logarithmic spiral structures in Nature, such as snail shells, 
galaxies, or eddies. Their iterative interactions reveal models for analogous interactions in 
Nature. We begin by defining the generating maps of an IFS Fractal, and go on to stating 
its existence and uniqueness. 

1.1 The Concept of IFS Fractals 

Definition 1.1 Let a 2- dimensional affine contractive mapping (briefly: contraction or con- 
traction map) T : C -> C be defined for all z G C as T(z) := p + (p(z — p) where p G C is the 
fixed point of T , and (p = \e d% G C is the factor of T , with A G (0, 1) the contraction factor 
of T , and 9 G (— 7r,7r] the rotation angle of T . 

In higher dimensions, an equivalent definition may be given using unitary rotation matrices 
R T R — I, R G IR nxn . Then the contraction maps take the form 

T(z) =p + XR(z -p), p,zeR n , AG (0, 1) 

This formulation shall be useful for our general bounding circle. 

Definition 1.2 A 2-dimensional affine contractive n-map iterated function system (briefly: 
IFS or n-map IFS) is defined as a finite set of contractions, and denoted as T := {Ti, . . . , T n }, 
n G N. We will denote the set of indices as Af := {l,...,n} ; the set of fixed points as 
V := {px, . . . ,p n }, and the set of factors as $ := {<pi, . . . , ip n }. 
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Definition 1.3 Let T = {T 1; . . . ,T n },n e N be an IFS. Define the Hutchinson operator H 
belonging to T as 

n 

H(S) : = |J T k (S), T k (S) := {T k (z) : z G S}, for any S C C 
fe=i 

Theorem 1.1 For any IFS with Hutchinson operator H, there exists a unique compact set 
FcC such that H(F) = F. Furthermore, for any compact Sq C C, the recursive iteration 
S n+ i := H(S n ) converges to F in the Hausdorff metric. 

Proof The proof follows from the Banach Fixed Point Theorem, after proving that H is 
contractive in the Hausdorff metric over compact sets [8]. □ 




Figure 4: Generation of a fractal by iterating a square (created by Scott Draves). 

Definition 1.4 Let the set F in the above theorem be called a fractal generated by an IFS 
with Hutchinson operator H (briefly: fractal, IFS fractal, attractor, or IFS attractor). Denote 
(Ti, . . . ,T n ) = (T) := F . We will say that F is a bifractal, trifractal, and polyfractal for the 
cases of n = 2,3 and any n > 3 respectively. In the special case of $k = 0, k 6 H we shall 
speak of Sierpinski fractals. 

Note that since by the above theorem, an IFS fractal is the limit set of a sequence that 
can also be considered a dynamical system, it is often referred to as an attractor of the IFS 
(or of H), which is a more general term. Mostly such attractors are self-similar and have 
non-integer (sometimes fractional) Hausdorff dimension, hence the term "fractal". There 
are cases when such an attractor is an ordinary Euclidean set. By introducing the last few 
terms, we wish to emphasize the "polygonal" nature of IFS fractals, especially apparent in 
the case of Sierpinski fractals. 
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1.2 Bounding Circles and Iterative Containment 

The relevance of a bounding circle lies in the fact that we may approximate the fractal 
via iteration of the Hutchinson operator over it, as illustrated by Figure 1 - this is the 
idea of "iterative containment" . Therefore a bounding circle serves as a tool for answering 
geometrical questions numerically to some e > accuracy. We hereby give an intuitive 
definition of a bounding circle. 

Definition 1.5 A circle C = (c, r) G C x IR + is called a bounding circle of the IFS fractal 
F, if it is contained in the disk determined by C , denoted as B(c,r) := {z G C : \z — c\ < r}. 

We shall refer to the iteration of sets numerous times throughout the text. This will mean 
the application of the Hutchinson operator H recursively to any set - primarily bounding 
circles - as in Theorem 1.1. The motivation being that for any B\ D F we have T k (F) C 
T k {B x ) (k G AT) and 

n n 

F = H(F) = |J T k (F) C |J T k {B x ) = H{B 1 ) =: B 2 

k=l k=l 

Since H is contractive in the Hausdorff metric over compact sets, then if B\ was compact, 
then so will B 2 be, and it will be closer to F in the metric. In practice, this will imply that 
taking any bounding set, such as a circle, we can improve it by one or several applications 
of H to new tighter bounding sets. 

Theorem 1.2 (Containment Theorem) If B C C is compact and H(B) C B then F C B. 
Proof The proof follows directly from Theorem 1.1. □ 

The above theorem provides an easily verifiable condition for whether a set bounds the 
fractal. This theorem is the key inspiration for our explicit bounding circle definitions. 

2 The General Bounding Circle for Polyfractals 

Let the radius function with respect to the IFS T = {Ti, . . . , T n } be denoted 

g(z) := max \p k — z\ (z G C) 

This is the radius of the smallest circle centered at z that contains all the fixed points p k G V . 
It can be shown that g is a continuous convex function, so it is a natural inquiry to look for 
its minimum - the center C* which gives the lowest g value. The minimizing circle (c*, f?(c*)) 
is called the minimum covering circle or minimal bounding circle. Determining it is the 
Smallest Circle Problem for a finite set of points in the plane, first investigated in modern 
times by Sylvester [15] in 1857. The minimum covering circle is known to be unique. Several 
algorithms exist for finding it, and the fastest run in linear time. We note the algorithms of 
Megiddo [12] and Welzl [17] specifically, both having linear runtime. 
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Theorem 2.1 For any c G C the circle 



(c,r(c))GCxM +) r(c):=^- 

A* := max |^|, /x* := max |1 — ^| 

zs a bounding circle of F. With the center c* = &rgtn.in zeC g(z) and radius r* := r(c*) this 
circle is minimal, and we will call this the general bounding circle of the IF S fractal F. 

Proof Considering the Containment Theorem, we see that (c, r(c)) must satisfy Tk(B(c, r(c))) C 
B(c,r(c)), k G N '. Let k be fixed and arbitrary. Then we can see geometrically that this 
containment is equivalent to 

\Tk{c) - c\ + X k r(c) < r(c) <=> |1 - <p k \ ■ |c - p k \ < (1 - A fc )r(c) 

Estimating the left-hand side by h*q(c) from above, the right-hand side by (1 — A*)r(c) from 
below, and the two estimates being equal, we have the desired inequality above. Therefore 
H(B(c,r(c))) C B(c,r(c)) so by the Containment Theorem we have that F C 5(c, r(c)), 
which also holds for (c*,r»). □ 

The above proof works in higher dimensional spaces for a corresponding IFS with maps 
of the form T k : lR n -> R n , T fc (2) = p fc + X k R k (z - p k ), RlR k = I and := ||J - A fe i2 fc ||, 
where || • || is the matrix norm induced by the Euclidean norm, or whatever norm was used in 
the definition of g (easy-to-compute compatible matrix norms with the vector norm are also 
sufficient). Thus in 3D we may talk about a general bounding sphere (c*,r*) of the above 
form, obtained from its corresponding minimal bounding sphere (c*, g(c*)). Algorithms do 
exist for finding the minimal bounding sphere of a finite set of points as well [9] . 

If we wish to obtain a good estimate of c*, we may work with the arithmetic or harmonically- 
weighted mean of the fixed points of the IFS. 

Ca ■= , c H - 



n J2 



k e{pk) 

For these special centers, we may find by the convexity of g that their value will be less 
than or equal to the arithmetic and harmonic mean of the g values of the fixed points, thus 
proving these centers reasonable. Calculating the g of both means, we may find whichever 
gives a lower value and use that center, rather than attempting the algorithmic calculation 
of the optimal c*. 

Philosophically speaking, we have reduced the problem of bounding an IFS fractal hav- 
ing an infinite number of points, to the bounding of a finite set of points, the fixed points 
of the IFS. Since we have pointed out specific efficient algorithms for finding c*, as well as 
reasonable explicit centers ca and ch, the definition of our general bounding circle can be 
considered an explicit one. 
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3 The Circumcircle for Trifractals 



The bounding of polygons by circles is in general formulaically non-explicit, since the min- 
imum covering circle of the vertices can only be determined algorithmically. When the 
polygon is a triangle however, such a circle - the circumcircle - may be determined explic- 
itly Since from a bounding perspective, Sierpinski polyfractals are essentially polygons of 
the fixed points of the IFS, their non-algorithmic explicit circular bounding seems to be a 
hopeless venture, for more than three maps. For the case of three or two maps however, 
such explicit circles may exist for any fractal. Since the Hutchinson of the circumcircle of a 
Sierpinski trifractal results in tangential iterates, it suggests a similar Apollonian definition 
for trifractals in general. This implies the following definition, illustrated by Figure 5. 

Definition 3.1 Let the circumcircle of a trifractal F = (Ti,T 2 ,T 3 ) be the smallest circle 
C = (c, r) G C x R + satisfying the conditions 

\T k (c) - c\ + X k r = r, A; = 1,2, 3 




Figure 5: A Sierpinski trifractal under rotational perturbations. 

Definition 3.2 Denote the dot and cross product of ' Zi i2 G C respectively by 

Z\ • z 2 ■= Re(zi)Re(z 2 ) + lm(zi)lm.(z2) , z\ x z 2 '■— Re(zi)Im(2; 2 ) — lm(zi)Re(z 2 ) 
Theorem 3.1 Let us introduce the following auxiliary variables 

«fc := h f> k = 1,2,3 

A := (c*3 - al)pi + (aj - a 2 3 )p 2 + (aj - al)p 3 

B := (\p 2 ? ~ N 2 )Pi + (N 2 - \Pi?)P2 + (hi 2 - \p 2 \ 2 )P3 
C := 2(p 2 -p l ) x (p 2 -p 3 ) 
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B A 
°° '~ C~\ r ° :— \ c ° ~ Pi 1 5 Cl :— 

£> := c • ci + ^(a 2 p 2 Xft + a 2 . p 3 x Pi + a 2 p x x p 2 ) 

T/ie circumcircle exists iffC^O, meaning iff Pi,2,3 non-collinear. 
If A 7^ ii /ias i/ie /orm 



c = c + cir 2 , r = t^— J-D - D 2 - \ci\ 2 rl 




If A = £/ien c = c and r = r§ja\. (Note: r can &e defined above with any pk-) 



Proof First we note that the circumcircle equations reduce to the conditions 



\p k - c| 2 - a 2 k r 2 = 0, A; = 1,2, 3 



Denoting pki : = Re(j9fc), Pfe2 := im (Pfc), ^ := Re(c), y := Im(c) and expanding these 
conditions, then multiplying each by the factors p%i — P21, Pu — P31, P21 — Pu respectively, 
and lastly summing the three equations, several terms drop out in the resulting equation, 
which we denote as E 1 = 0. We may do similarly with the factors P32—P22, P12—P32, P22—P12, 
resulting in an analogous equation E 2 = 0. Taking E\ + E 2 i = and collecting terms, we 
get that Ar 2 + B — Ci ■ c = 0, which gives the formula for c. To find r, we expand the 
circumcircle condition say for k — 1, resulting in 



When ci = 0, we get the stated formula for r. When c\ 7^ 0, we may solve this equation for 
r 2 . Upon some calculation, we deduce that the coefficient of the r 2 term in the equation is 
2D. The quadratic formula gives the desired expression for r, by choosing the negative sign in 
±, since we have defined the circumcircle to be the smallest circle satisfying the conditions. □ 

We remark that when = for some k, then ak = 1, so by the corresponding circumcircle 
condition, we get that \pt — c\—r. This means that when dk = 0, the fixed point pk lies on 
the circumcircle. For Sierpinski trifractals, all three fixed points lie on the circumcircle, as 
expected. 



Similarly to the case of three maps, we will define the two-map circumcircle in an Apollonian 
manner, to "optimally" satisfy the condition of the Containment Theorem. We derive it in 
quite a different manner however. Having any bounding circle to the fractal, we may get 
a potentially smaller one, by applying the Hutchinson operator to it, and then taking the 




4 The Circumcircle for Bifractals 
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outer tangential circle as in Figure 6 below - this will be calculated using the map M. So 
the two-map circumcircle will then be defined as the fixed point - or fixed circle - of M, 
containing the fractal by its definition and the Containment Theorem. 

If (a, r>) and (b, re) represent the circles received from a bounding circle under the ac- 
tion of Ti and T 2 respectively, then finding their outer tangential circle (c, rc) will usually 
result in a tighter bounding circle than the original. 




Figure 6: The Circumcircle for Bifractals. 



From Figure 6, the following relations can be seen 

c = ^(( a _ t aUab) + (p + t b uab)), r c = \{r A + r B + \b - a\), Uab 



2 yy « ■ \ - — //> - 2 V A 13 1 Iy ' ™ \b-a\ 

These relations suggest the following transformation. Given a circle C = (c, r) in the plane, 
let M :— (mi,m 2 ) :CxI->Cxl be defined as 

. T 1 (c)+T 2 (c) \ 2 -\iT 2 (c)-Tx(c) 
mi(c, r) := h r 



T 2 c -T! c 



777.2(0, rj := rH 



Theorem 4.1 Lei i/je circumcircle of a bifractal be defined as the unique fixed point of M , 
which is given by the following (c, r) 6 C x 1R + 

_ _ (1 - f)(l - ¥>0P1 + (1 + f)(l - ¥? 2 )p 2 

|P2 -Pl| 



(1_ I/ )(1_ ¥ , 1 ) + (1 + I/ )(1_^ 2 ) 
|1 - <^i| • |1 - (f 2 \ 



(l-A)-|(l-i/)(l-^) + (l + i/)(l-^ 2 )| 
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Ai + A2 A2 — Ai 

: " 2 ' V := 2(1 -A) 



Proof The fixed point equation M(c, r) = (c, r) implies that rri2(c,r) = r which implies 
that r = ||T 2 (c) — Ti(c)|/(1 — A). By plugging this into the equation mi(c,r) = c and some 
further manipulation, we get the above formula for c. To get the formula for r, according to 
the above equation it is sufficient to calculate T 2 (c) — Ti(c) in terms of pi^ and <^i )2 . Since 
we have calculated c, this is a matter of some further algebraic manipulation. The derivation 
shows that the fixed point of M not only exists, but that it is also unique. We note that c 
is the complex combination of p\ and p2- □ 

5 Remarks 

We hereby wish to make a brief but important note. One need not be satisfied working with 
any given bounding circle C = B(c,r) of F, since it can be easily improved. Having such a 
circle and taking some finite L-level iterate of it by the Hutchinson operator, ie. H L (C), we 
can take the smallest enclosing circle (c', r') of the n L centers H L ({c}), and then (c', r' + X^r) 
will be a tighter bounding circle of F. For large enough L, we can get within any s > 
accuracy of the fractal. So we do not need any complex algorithms to arrive at the tightest 
bounding circle - the minimum bounding circle and sphere algorithms resolve the problem. 

Upon considering the various circles, one may contemplate whether there is any good reason 
to introduce these circumcircles for the case of two or three maps. There is however, in that 
they may be potentially tighter than the general bounding circle. According to our numerical 
experiments, in most cases the circumcircles are somewhat tighter, yet the general bounding 
circle still remains competent. In a few cases, it is even tighter than a circumcircle, so for 
applications, we advise the programming of all three circles. 




Figure 7: The circumcircle (red) vs. the general bounding circle (blue). 



Throughout our paper, we have emphasized the potential of bounding circles to resolve 
virtually any problem in the plane numerically, thus being a "panacea" . We conclude with an 
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illustration of the numerical solution to the fractal-line intersection problem, which suggests 
analogous solutions to other geometrical questions. In this and other specific cases, the 
algorithms may be designed ideally for the lowest possible algorithmic complexity. 




Figure 8: Intersecting a fractal with a line numerically. 



6 Conclusion 

We have introduced three kinds of explicitly-defined bounding circles for IFS fractals with 
similitude contractions. The general bounding circle encloses polyfractals, and it even gen- 
eralizes to higher dimensions. The circumcircle of trifractals generalizes the analogous circle 
for triangles, hinting at a novel polygonal view of IFS fractals. The circumcircle of bifractals 
generalizes the concept of the interval for Cantor sets on the real line. Indeed tri- and bifrac- 
tals may hold a special place in the theory of Fractal Geometry, and may require as much 
attention as triangles in Euclidean Geometry. Cyclic polyfractals, where the fixed points lie 
on a circle, may be analogously rich in research potential. 

Our investigations indeed intend to spawn future mathematical research, and our paper has 
several inspirational half-thoughts that may be elaborated upon. To point out some, the 
algorithm for the fractal-line intersection problem may have an optimal design, whether 
in 2D or 3D. There is a wide variety of other geometrical problems, having an optimal 
algorithmic solution using bounding circles. On a more theoretical level, cyclic polyfractals 
may also possess an explicit bounding circle, expressible in terms of the circle along which 
the fixed points lie. The ultimate question of the explicit tightest bounding circle of any IFS 
fractal lingers. Extending results to higher dimensions remains a continual challenge. 

We hereby express our gratitude for the support of the Faculty of Mathematics at the Uni- 
versity of Waterloo, and primarily that of E.R. Vrscay who is supported by NSERC. 



12 



References 



[1] M. F. Barnsley. Fractals Everywhere. Academic Press, Boston, second edition, 1993. 

[2] M. F. Barnsley and D. P. Hardin. A mandelbrot set whose boundary is piecewise smooth. 
Transactions of the American Mathematical Society, 315(2) :641-659, 1989. 

[3] D. Canright. Estimating the spatial extent of the attractors of iterated function systems. 
Computers and Graphics, 18(2):231-238, 1994. 

[4] S. Dubuc and R. Hamzaoui. On the diameter of the attractor of an ifs. Technical report, 
C.R. MATH. REP. SCI. CANADA, 1994. 

[5] A. Edalat, D. W. Sharp, and R. L. While. Bounding the attractor of an ifs. Technical 
report, Imperial College, 1996. 

[6] K. J. Falconer. Fractal Geometry. Wiley, Chichester, second edition, 2003. 

[7] J. C. Hart and T. A. DeFanti. Efficient antialiased rendering of 3-d linear fractals, 
volume 25 of Computer Graphics, pages 91-100, Las Vegas, 1991. SIGGRAPH. 

[8] J. E. Hutchinson. Fractals and self similarity. Indiana University Mathematics Journal, 
30:713-747, 1981. 

[9] T. Larsson. Fast and tight fitting bounding spheres. In Proceedings of The Annual 
SIGRAD Conference, pages 27-30, November 2008. 

[10] B. B. Mandelbrot. The Fractal Geometry of Nature. Freeman, San Francisco, 1982. 

[11] T. Martyn. Tight bounding ball for affine ifs attractor. Computers and Graphics, 
27(4):535-552, 2003. 

[12] N. Megiddo. Linear-time algorithms for linear programming in r3 and related problems. 
SIAM Journal on Computing, 12:759-776, 1983. 

[13] J. Rice. Spatial bounding of self-affine iterated function system attractor sets, pages 
107-115. Graphics Interface, 1996. 

[14] D. W. Sharp and R. L. While. A tighter bound on the area occupied by a fractal image. 
Technical report, Imperial College, 1999. 

[15] J. J. Sylvester. A question in the geometry of situation. Quarterly Journal of Mathe- 
maitcs, 1:79, 1857. 

[16] E. R. Vrscay. Mandelbrot sets for pairs of affine transformations in the plane. J. Phys. 
A: Math. Gen., 19:1985-2001, 1986. 

[17] E. Welzl. Smallest enclosing disks (balls and ellipsoids). New Results and New Trends 
in Computer Science, Lecture Notes in Computer Science, 555:359-337, 1991. 



13 



